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A semiclassical theory of single and multi-mode lasing is derived for open complex or random 
media using a self-consistent linear response formulation. Unlike standard approaches which use 
closed cavity solutions to describe the lasing modes, we introduce an appropriate discrete basis of 
functions which describe also the intensity and angular emission pattern outside the cavity. This 
constant flux(CF) basis is dictated by the Green function which arises when formulating the steady 
state Maxwell-Bloch equations as a self-consistent linear response problem. This basis is similar 
to the quasi-bound state basis which is familiar in resonator theory and it obeys biorthogonality 
relations with a set of dual functions. Within a single-pole approximation for the Green function 
the lasing modes are proportional to these CF states and their intensities and lasing frequencies 
are determined by a set of non-linear equations. When a near threshold approximation is made to 
these equations a generalized version of the Haken-Sauermann equations for multi-mode lasing is 
obtained, appropriate for open cavities. Illustrative results from these equations are given for single 
and few mode lasing states, for the case of dielectric cavity lasers. The standard near threshold 
approximation is found to be unreliable. Applications to wave-chaotic cavities and random lasers 
are discussed. 
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I. INTRODUCTION 

A long-standing problem in laser theory is the formula- 
tion of a model for lasing which correctly treats the open- 
ness of the lasing medium/cavity and the non-linearity of 
the coupled matter- field equations. This mathematical 
challenge has become of great relevance with the current 
high interest in complex or random lasers, for which the 
mode geometry is not fixed by the placement and orienta- 
tion of mirrors and one typically has multi-mode lasing 
behavior. In this case many spatially complex internal 
modes contribute to the external emission pattern and 
we currently lack any theory to predict the directionality 
of emission patterns and to understand and predict the 
output power as a function of pump strength. One class 
of systems of particular interest in this regard arc dielec- 
tric cavity lasers with complex and sometimes chaotic ray 
dynamics 111 ; another class is the so-called "random" 
lasers 0, U in which light undergoes diffusive motion 
within the gain medium. In the former case the laser will 
typically have relatively high-Q modes whose lasing prop- 
erties are nonetheless determined by complicated compe- 
tition between spatially complex modes. The latter part 
of this work will focus on the dielectric cavity case al- 
though the formalism we are developing should be useful 
for both random lasers and for certain conventional lasers 
in which mode competition is important. We believe the 
current approach solves the problem of treating the open- 
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ness of the cavity in the simplest possible manner, and 
makes clear the connection between the resonances (or 
quasi-bound states) of the passive (cold) cavity and the 
lasing modes of the active cavity. In particular if a pe- 
riodic or multi-periodic solution of the Maxwell-Bloch 
(MB) equations exists, then the lasing modes and fre- 
quencies are determined by a set of self-consistent inte- 
gral equations the kernel of which is the Green function 
of the linear problem with outgoing wave boundary con- 
ditions. When the modes are relatively high-Q a simple 
approximation to this Green function implies that the 
spatial modes are given by a certain set of states which 
we refer to as constant flux (CF) states. The effects of 
spatial-hole burning are then described by interactions 
between these modes, which unlike previous theories, are 
modes which exist in all space and can thus be used to 
predict output power and emission patterns from com- 
plex or random two and three-dimensional laser cavities. 
Throughout this work we employ the semiclassical de- 
scription of lasing implied by the MB equations, so the 
effects of field quantization, such as quantum fluctua- 
tions, are not included. 

The modal description of lasing assumes that the laser 
is in a regime in which the steady-state electric field in 
the pumped medium E(x, t) (which we will henceforth 
refer to as the "cavity") has a finite number of frequen- 
cies and hence by Fourier transform can be seen as a sum 
of non-linear modes E tl (x,n fl ). It is now well known 
that lasers can have chaotic temporal dynamics which 
cannot be described by a finite number of frequencies 
000 For clarity we wish to emphasize that wc 
are not treating lasers which exhibit dynamical chaos 
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in this sense; the "chaotic lasers" we are treating are 
those with complex modal patterns which can be related 
to the chaotic motion of light rays in the geometric op- 
tics limit. A term for this type of laser consistent with 
usage in the field of quantum chaos is "wave-chaotic". 
The temporal dynamics we are treating is conventional 
multi-mode (MM) lasing, typical of the vast majority of 
lasers. There have been a number of interesting exper- 
iments for which the emission patterns from deformed 
(non-spherical or cylindrical) dielectric cavity lasers have 
not been explainable in terms of conventional "whisper- 
ing gallery" modes, but instead required careful analysis 
of different modal patterns with different ray-optical in- 
terpretations jlGlillllEEQEIIii- Furthermore 
several experiments have found a dramatic variation of 
the output power for a given pump strength with the 
shape of the laser cavity j^, Hol | . Our current formalism 
is proposed to describe and predict the results of such 
experiments, which are inadequately treated by standard 
approaches. 

There are several methods to predict or describe lasing 
modes and emission in this case of multi-periodic steady 
state solutions. The simplest approach is to work with 
the "cold cavity" (CC); i.e. to find the electromagnetic 
solutions for the Helmholtz equation describing light in 
the cavity neglecting gain and non-linear effects. If the 
lasing modes have high Q values then the non-linear so- 
lutions are often very similar spatially to these CC modes 
[iH IT^ and from the experimentally observed MM las- 
ing frequency spectrum one can associate the lasing state 
with a sum of cold cavity modes (or a single mode if only 
one frequency is present). Within this CC approximation 
a further approximation is to treat the cavity as closed 
and work with hermitian modes of a perfectly confining 
cavity. The Gaussian modes of Fabry-Perot resonators 
shown in textbooks are solutions of this type. 

Another more common approach within the CC ap- 
proximation is to approximate the lasing modes by 
the quasi-bound (QB) states or resonances of the cavity. 
These states arc defined as the solutions of the linear 
Maxwell wave equation (in the absence of gain) which 
satisfy the boundary conditions at the cavity boundary 
and only have outgoing waves at infinity. This outgoing 
wave boundary condition can only be satisfied for discrete 
complex values of k = q — 17 (with 7 > 0) which means 
that the outgoing spherical waves will grow as exp[7r] 
as r — ► 00. Such modes are not hermitian and are not 
orthogonal to one another (though a modified orthogo- 
nality relation can be defined in separable non-hermitian 
problems 0, l2p|). The quantity 7 determines the Q- 
value of the resonance which can then be used in formu- 
lating the laser theory. The Fox-Li method [2l| dating 
from the early days of laser theory is one technique for 
finding such quasi-bound modes; there are a number of 
modern methods for doing this as well [1 Ej I2I l24l . 
This theoretical approach has been widely applied to an- 
alyze post facto the modes of wave-chaotic or random 
lasers (see pj for references). Its obvious drawback is 



that the method is not predictive. First, a larger set 
of QB modes is found and a subset, the lasing modes, 
are chosen as those which "look" like the experimental 
results. Furthermore, there is no means within this ap- 
proach to understand and predict the output power of 
the laser. 

The standard approach to go beyond the CC approxi- 
mation is to use the simplest semiclassical lasing theory, 
describing a pumped medium of two-level atoms coupled 
to light within a cavity, known as the Maxwell-Bloch 
(MB) equations. Our theory below is based on analy- 
sis of these MB equations. The MB equations can be 
analyzed to find single-frequency solutions, usually ap- 
proximated by single modes of the cold cavity. One then 
finds that it is the highest-Q cold cavity mode which 
lases first and that its frequency is shifted from the cold 
cavity frequency towards the atomic transition frequency 
but not by a large amount if the Q of the cavity is suffi- 
ciently high 0, 0, 0, [2(| . The spatial field distribution 
is typically assumed unchanged from the CC mode which 
can be approximated by the closed cavity mode or the 
corresponding CC resonance. The complication in the 
theory comes when one attempts to describe multi-mode 
lasing. First the general non-linear MB equations cannot 
be solved analytically using a modal expansion of the 
electric field. Exact numerical solution of the equations 
can be done for some cases and is useful ; however this 
becomes computationally intractable in the short wave- 
length limit of interest here and it is difficult to extract 
qualitative physical ideas from such an approach. A nice 
method due to Haken dates back to the early days of 
laser theory; one limit of our theory is an extension of this 
approach. The Haken method expands the electric field 
which solves the MB equations as a sum of CC modes 
and writes an equation of motion for the amplitude of 
each mode which can be reduced to a (constrained) lin- 
ear equation for the modal intensities in the near thresh- 
old approximation. The major drawback of the Haken 
approach is that it is formulated in terms of the modes 
of the ideal closed cavity and calculates only the inter- 
nal field intensities of these modes; thus without some 
ad hoc assumptions it doesn't predict output power or 
directional emission patterns. This is done to exploit the 
orthogonality of closed cavity modes which simplifies the 
mathematical description. For comparison to our results 
below, the Haken near threshold MM lasing equations for 
the steady state electric field intensities /„ are 

where k^, Do, g{Qv) are the cavity decay rate, pump 
strength, and laser gain profile respectively (in appro- 
priately scaled units) and 

= /Vss|v? /1 (x)| ! V(x)| 2 (2) 

is a matrix which describes interactions between CC 
modes ip^(x), tp v (x). The complexity of modal solutions 
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in space thus determines the lasing amplitudes through 
the properties of this matrix. This equation was derived 
by Haken and Sauermann in 1963 to describe the effects 
of spatial hole-burning: the fact that lasing modes de- 
plete the inversion in a spatially-varying manner which 
can then allow many modes to lase in steady-state when 
the pump is sufficiently strong. A statistical analysis of 
certain aspects of these equations for the case of closed 
cavities described by random matrix theory was per- 
formed by Misirpashaev and Beenakker some time ago 

a. 

Note that while Eq.QJappears to be a simple inhomoge- 
neous linear equation for the {I u }, it cannot be trivially 
inverted to yield these intensities since we have the ad- 
ditional constraint that I u > 0,V„. The equation thus 
needs to be solved by an iterative procedure. In a recent 
work [2^ we describe this procedure and used it along 
with some ad hoc assumptions to calculate multi-mode 
output intensities for deformed dielectric cylinder lasers, 
making contact with the experimental results of Ref . . 

In the current work we present a formalism which gives 
a generalization of Eq. ^ and which avoids these ad hoc 
assumptions. The basic innovation of the approach is 
to formulate the solution of the MB equations as a self- 
consistent linear response problem and write this solution 
in terms of the Green function of the wave equation with 
the boundary condition of constant outgoing flux at in- 
finity. Using this Green function the lasing solutions and 
frequencies are determined by an integral equation which 
can be used to describe multimode lasing in both high-Q 
and low-Q cavities, and in random media (acting as a cav- 
ity). To describe complex cavities with high-Q modes it is 
useful to represent this Green function in terms of a new 
set of linear eigenfunctions which satisfy the wave equa- 
tion with real wavevector but purely outgoing boundary 
conditions; we refer to these functions as constant flux 
(CF) states. They are similar to the resonant solutions 
of the cold cavity in that they have a complex wavevec- 
tor inside the lasing medium, but they differ in that they 
have a real wavevector outside the medium and hence 
give a well-defined field at infinity. Like the quasi-bound 
states, the CF states are not orthogonal, and the Green 
function involves both the CF states and their biorthogo- 
nal partners which represent states of constant incoming 
flux from infinity. The use of biorthogonal pairs of res- 
onator states is well established in resonator theory, and 
it well is known that the two states correspond to dif- 
ferent directions of propagation through the apparatus 
(see Ref. H3], PP 847-857). Our CF states are just a 
generalization of this idea. The new feature here is that 
a simple "single-pole" approximation to the Green func- 
tion implies that the lasing modes are proportional to 
the CF states and leads to a generalization of the Haken- 
Saucrmann multi-mode lasing equations in terms of CF 
states. Iterative solutions of these equations predict the 
multi-mode lasing states both inside and outside the cav- 
ity and hence provide a predictive theory of output power 
and directional emission from complex cavities. 



The paper is organized as follows. In Section [H] we 
derive the self-consistent integral equation for the lasing 
modes assuming that a multi-periodic solutions exists. 
To do this we introduce a Green function for the inhomo- 
geneous wave equation within the cavity. In Section IlIII 
we discuss the boundary conditions on this Green func- 
tion and its spectral representation. We show that in or- 
der to satisfy the outgoing wave boundary conditions this 
spectral representation is expressed in terms of two sets of 
biorthogonal functions (the CF states). The approxima- 
tion corresponding to standard multi-mode lasing theory 
is to approximate the Green function by the contribu- 
tion of the single-pole nearest to the lasing frequency. 
This implies that the lasing modes are proportional to a 
single CF state. We discuss examples of CF states for 
a slab (ID), cylinder (2D) and deformed cylindrical di- 
electric cavity. In Section IIVI we derive the multi-mode 
lasing equations which follow from this approximation, 
leading to linear and non-linear generalizations of the 
Haken-Sauermann equations, valid for the open cavity. In 
Section [V] we discuss iterative methods for solving these 
equations for the modal intensities and the lasing fre- 
quencies and then present solutions for the single and two 
mode cases, beyond the standard Haken-Sauermann near 
threshold approximations. Our results suggest that the 
near threshold approximation introduces significant error 
and greatly overestimates the number of lasing modes at 
a given pump power. We summarize our results in Sec- 
tion I VII and provide further detail on CF states in the 
Appendices. 



II. DERIVATION OF THE SELF-CONSISTENT 
SEMICLASSICAL LASER EQUATIONS 

Following the standard semi-classical laser theory 0, 
W\ . we start with the Maxwell-Bloch equations (MB) in 
the form 



c c 



P+ + (PH 



ih 



ED 



(3) 



(4) 



D 



^(D -D)--E[(P+r-P+] (5) 



This is a set of non-linearly coupled spatio-temporal 
partial differential equations for the electric field am- 
plitude E(x,i), the macroscopic polarization P(x,t) = 
n a(gpi2(x, t) +g*p2i(x, t)) = P + + P~ and the inversion 
D(x, t) = n a (p 22 (x 1 t) - Pll (x, t)). Here, P+ = (P")* = 
n a gpi2- The parameters entering the equations are as fol- 
lows: Dq is the external pump strength, g is the dipole 
moment matrix element, 7j_ and 7|| are phenomenological 
damping constants for the polarization and the inversion, 
respectively, and c is the speed of light. Note that E, D 
and P are real valued fields. 

We will make the following assumptions: 
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• We focus here on a scalar field E defined in two 
space-dimensions x = (x,y). For instance, in the 
case of a dielectric cylinder laser with arbitrary 
cross-section, E will denote the z-component of the 
electric or magnetic field for k z = modes (see 
Ref. pj for the justification of this model in the 
case of a dielectric cylinder laser). 

• We assume a uniform, homogeneously broadened 
atomic medium with density of atoms n a , atomic 
transition frequency tu a and the quantum mechan- 
ical density matrix pij(x,t). 

We decompose P = Pl + Pnl, into a linear compo- 
nent Pl (for instance the non-resonant response of the 
substrate material), and a non-linear resonant compo- 
nent Pnl due to the gain medium, here taken to be the 
uniformly distributed set of two-level atoms. Let 

P L (x,t) = J dt'x(x,t-t')E(x,t') (6) 

Decomposing the fields P = Pl + Pnl and E as 

F(x,t) = F + (x,t)e- iUat +F-{x,t)c tu " t (7) 

where F = {E , Pl, Pnl) and using the slowly varying 
envelope approximation and the rotating wave approxi- 
mation, Eq. |3 becomes 



2iuj n 



V 2 e+^(e + 47rp L ) + ^-(e + 47rp L 



-PNL (8) 



The envelope equations corresponding to Eqns. |2E1 be- 
come 



d t (n 2 e) 

Pnl 
D 



— [u? a h 2 + c 2 V 2 ] e + 2iTTiv a pNL (9) 

ZbJ a 



-J±Pnl + -zbD 
in 



(10) 



7|| (A, -D)--^ (Wnl ~ Pnlc*) (11) 



Here, p NL (x,t) = P+ L (x,t), p L (x,t) = P^{x,t), 
e(x,t) = E + (x,t), and n 2 stands for the convolution 
operator J dt'n 2 (x,t — t'). The time-dependent index 
of refraction n(x,t) is related to x( x ,t) by n 2 {x,t) = 
1 + 4nx(x,t). In the derivation of these equations we 
have assumed that the residual time-dependence of the 
envelopes is much slower than cu a and that u a 3> 7_l and 
dropped the terms proportional to e,pNL,^aPNL,PL- 

The properties of the index of refraction n(x, t) of the 
background medium will depend on the particular com- 
plex lasing system under investigation; for dielectric cav- 
ity lasers it will be constant in the cavity and unity out- 
side, for a random laser it will vary randomly in space 
within some region and then fall off to unity at the edges 
of the medium. For the wave-chaotic dielectric cavity 
the "randomness" comes from the scattering at an irreg- 
ularly shaped boundary. We will also replace Pnl by p 
henceforth. 



We assume now that the fields p(x,t) and e(x,t) are 
multi-periodic in time 



p(x,t) 



Pf,{x)e 



—iSl^t 



e(x,t) = 5^tf M (sc)e- 



(12) 

in the steady state. In contrast to previous approaches 
IE HE HH HE HE IS- we leave the spatial functions p^(x) 
and *5> ^(x) to be unknown functions and to be the 
unknown lasing frequencies to be determined. Such a 
solution with a finite number of discrete frequencies (a 
"multi-mode" lasing solution) is only possible if the in- 
version is approximately stationary |35j . Henceforth, we 
will assume this to be the case. From Eq. (|10|) . the ansatz 
Eq. (|12|) and the time-independence of the inversion we 
obtain 

a 2 Dfx) 

M*) = ^T .-o , .. ( 13 ) 



and 



D(x) 



ih —iQfi + 7_i_ 



i + ^£ M s(^)l*>)P 



(14) 



Here, g(il) = 7^/(f2 2 + 7j_) is the laser linewidth and 
e c = h y / , y± r y\\/2g gives the typical electric field scale. 
Henceforth we will measure the lasing modes and polar- 
ization, \& ^(x),p^(x) in units of e c . We can then rewrite 
Eq. in the form 

[2iu a d t n 2 + u 2 a n 2 + c 2 V 2 ] e(x, t) = -4irtulp(x, t) (15) 

Substituting the results for p^x) into Eq. I|15fl . the 
pumped atomic medium becomes equivalent to a multi- 
periodic forcing term on the right side of the Eq. i|15fl . 
given by 



-47TU Z aP (x, t) = -47TW 2 



29 



Do 



ifil + E,5(^)l*M(^)l 2 

■* M (x)e-* n "* 



*-f + 71. 



(16) 



where the polarization which is the source for e(x, t) = 
^2 ^ ll {x)e~ in ^ t is itself a non-linear function of {^j. 
Thus Eq. (|T3|) gives a self-consistent set of equations for 
the mode amplitudes ^^(x) and lasing frequencies f2 M . 

To derive from Eq. I|15|) an integral equation for each 
^^(x), we first Laplace transform Eq. i|15|) to the fre- 
quency domain: 

(oj a + Lu) 2 n 2 (uj)e(x, u) + c 2 S7 2 e(x, uo) = —^■klo 2 p(x, uo) 

(17) 

Here, n 2 (uj) = 1+47tx(cj), the Laplace transforms are de- 
fined by p(co) = ^ J °° dte^+^pit), and we only retain 
the infinitesimal imaginary part of the frequency, e, when 
it is needed for convergence. Note that the real frequency 
variable u> corresponds to the residual time-dependence 



5 



of p(x,t) after removing the rapid variation at frequency 
uj a 3> uj, and hence we have used uj 2 + 2uiuj a « (u> a + u>) 2 . 

Initially we treat the right hand side of Eq. (|17J) as 
a given source and introduce a Green function for the 
inhomogeneous equation, which ensures that the solution 
e(x, uj) satisfies the appropriate boundary conditions, i.e. 
that the field has only outgoing contributions. Assuming 
for the moment that such a Green function exists we can 
immediately write a formal solution to Eq. fll 51) within 
the cavity 



e(x, uj) 



2iri(uj — il M + ie) 



dx' G(x,x'\uj)p(x',uj) (18) 



cavity 



The integral on the rhs of Eq. I|18(l is only over the cavity 
because the source (pumped medium) is zero outside of 
the cavity. The cavity Green function satisfies the equa- 
tion 

[V 2 + + uj a ) 2 }G(x, x'\ui) = S 3 (x - x') (19) 

which formally inverts Eq. i|17|) to yield Eq. i|18|) for the 
electric field component e(x,u) for x within the cavity. 
Outside the cavity the wave equation changes and a dif- 
ferent Green function would be needed to invert the equa- 
tion. However as the cavity Green function determines 
the solution on the boundary of the cavity, one can just 
use the outgoing wave boundary condition to fix the so- 
lution outside. This will be done through our definition 
of the CF states below. 

Both sides of Eq. <|18|) involve sums over the modes p 
weighted by the factor \/2iri(uj — + ie) which will give 
the multi-periodic time dependence of the solution when 
it is inverse Laplace transformed. The analytic structure 
of the relevant Green function can be shown to produce 
no additional harmonic dependence in time as t — > oo 
so that it is possible to equate the residues of these real 
poles to yield an integral equation for each of the modes 



47tuj 2 a g 2 Do 



f G(x,x'\n^,(x') 

J cavity 



(20) 



This set of non-linear integral equations for the lasing 
modes is completely general assuming a multi-periodic 
solutions exists. If the modes of the cavity are sufficiently 
high-Q that the imaginary part of their frequency is much 
smaller than the typical mode spacing then one can intro- 
duce an approximation for the Green function of Eq. I|20|) 
which leads to generalizations of the conventional multi- 
mode lasing equations (see below). However, it is possible 
to work directly with Eq. (|20|) when, as in a random "dif- 
fusive" cavity, the assumption of well-separated high-Q 



modes is not correct. In this case different approxima- 
tions to the Green function, similar to those made for 
disordered electronic systems, would be more appropri- 
ate and would yield composite lasing modes generated 
from many very broad cold cavity resonances. We will 
not pursue this direction in the current work, but wish 
to point out the feasibility of such an approach. 



III. BIORTHOGONAL CONSTANT FLUX 
STATES 

A. General Definition and Properties 

To express the lasing solutions ^^(x) as an expansion 
in a set of linear cavity modes we now introduce a spectral 
representation of the Green function G(x, x'\ui) of the 
form: 



G{x,x'\lu) = 



(p m (x,uj)ip* m (x',uj) 
n 2 [{k + k a ) 2 - (k m + k a ) 



(21) 



where k = uj/c,k a = uj a /c and the functions in the 
numerator are chosen so that G satisfies the outgoing 
wave boundary conditions. It is here that the non- 
hcrmitian nature of these boundary conditions requires 
a change from the standard approach utilizing closed 
cavity (hcrmitian) modes. We needed to introduce in 
Eq. (|21|) two sets of functions which are biorthogcmal 
{ip m (x, uj)}, {(p m {x, uj)}. The set {ip m (x, uj)} satisfies the 
eigenvalue equation 

-V 2 ^ m (x,w) = n 2 (k m + k a ) 2 ip m (x,uj) (22) 

defined in the cavity, T>. While the Green function of 
Eq. (|21|l is only used to find the electric field inside the 
cavity (and would give an incorrect answer for the electric 
field outside), one can define the functions ip m outside 
the cavity as well. Outside the cavity they satisfy the 
free wave equation with the fixed external wavevector, 
k a + k, (not k a + k m ): 



-V 2 (p m (x,uj) = (k a + k) 2 ip m (x,uj), 



(23) 



with the outgoing wave boundary conditions expressed 
in d-dimensions by: 



r^-V I 2 <p m (r -> oo, uj) ~ e ^ +k > , 



(24) 



equivalently it can be expressed linear homoge- 

neous boundary condition on the function $ m (a;) = 
r^ d ^ 1 ^ 2 ip m (r) in the form din /dr = i(k a +k) as r — > oo. 

Thus each ip m satisfies a different differential equation 
inside and outside the cavity but are connected by the 
continuity conditions 

V>m\dT>- = l Pm\dV+ j d n (p m \ dT ,- = d n ip m \g v + (25) 

at the boundary of the cavity (here &D denotes the 
boundary of T> and d n is the normal derivative on dT>). 
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This choice guarantees that in the approximation devel- 
oped below, in which each lasing mode is proportional 
to a single ip m , the electric field will also be continuous 
at the dielectric boundary, satisfy the free wave equation 
outside, and have only outgoing components at infinity. 
The spectral representation, Eq. Ij21|l is only meaningful 
for x, x' within the cavity and on its boundary. We will 
focus on the typical case in which one can impose the 
outgoing wave boundary condition directly on the cav- 
ity boundary (as in the slab and cylinder examples given 
below). In this case the eigenvalues (k m + k a ) 2 are deter- 
mined by linear homogeneous boundary conditions on a 
finite region, and are complex due to the complexity of 
the logarithmic derivative at the boundary. 

Note however, that this condition differs subtly from 
the condition defining the quasi-bound (QB) states, for 
which the complex eigenvalue k m itself would appear in 
the logarithmic derivative at the boundary, and hence in 
the external outgoing wave. As a result, as noted earlier, 
the QB states grow exponentially at infinity and carry 
infinite flux outwards. The states tp m we have just de- 
fined have real wavevector at infinity and carry constant 
flux; hence we refer to them as constant flux (CF) states. 

The CF states have appeared naturally from the con- 
dition that the Green function of Eq. (|18[) satisfy the cor- 
rect radiation boundary condition; this allows us to for- 
mulate a lasing theory valid for an open cavity in terms 
of these states. In the simplest approximation, to be de- 
veloped below, the lasing modes are just single CF states 
(those above threshold for a given pump power), each 
scaled by an overall intensity factor determined by solv- 
ing Eq. (EOJ, or its near-threshold approximation, which 
turns out to be a generalization of the Haken-Sauermann 
equations. Note that CF states have complex wavevec- 
tor inside the cavity, and are amplified while traveling in 
the gain medium, but real wavevector and conserved flux 
outside as one expects for lasing modes (see Fig|21below). 
We show in the Appendix that for high-Q modes the CF 
wavevector inside the medium is very close to that of the 
corresponding QB state, justifying the use of QB states 
to calculate the Q value of lasing modes, as is often done. 

As mentioned already, due to the non-hermitian nature 
of the associated eigenvalue problem we need to introduce 
a second set of functions, which satisfy different boundary 
conditions, in order to represent the Green function of 
Eq. J2H- The functions {(p m (x,uj)} satisfy 

~V 2 (p m {x,uj) = n*{uj) 2 (k m + k a ) 2 (p m {x,uj) (26) 

with the continuity conditions l|25|) but with the incoming 
wave boundary condition 



they satisfy biorthogonality with respect to the following 
inner product 



,(d-i)/2 



-i{k a +k)r 



(27) 



It can easily be shown |3y| that = k 2 ^ but that 
(p m (x,u>) 7^ Lp* m in general. (We will see in an exam- 
ple below that when (p m {x, u) is real for real frequencies, 
then (p m (x,oj) does equal tp^)- 

The key property which motivates the use of the two 
sets of functions in the spectral representation is that 



(((p m \fn)) 



d 2 : 



-,(p* m (x,u)(p n {x,w) = Tjm^Sj, 



(28) 

The normalization factor r/ rn can be set equal to unity 
for any specific choice of the frequency w, but it is useful 
to retain it explicitly in some contexts to denote that the 
normalization depends on uj. Note that we have here the 
usual hcrmitian inner product, except that it is defined 
between eigenfunctions living in adjoint spaces. These 
pairs of functions with eigenvalues related by complex 
conjugation are called biorthogonal pairs or partners. 

The boundary condition (|27|l clearly has the meaning 
that the adjoint CF states have constant incoming flux at 
infinity. For these states one can view the source as being 
at infinity and emitting radiation which impinges on the 
dielectric medium and then decays inwards at the same 
rate as the CF states grow while moving outwards from 
the origin. In other words the dielectric medium acts as 
an amplifying medium for outgoing CF states and as an 
absorbing medium for ingoing (adjoint) CF states. 



B. CF states and multi-mode lasing 

The importance of CF states becomes clear when we 
consider the self-consistent integral equations for the las- 
ing modes (f^Uf) using the CF spectral representation of 
the Green function. We can write Eq. (|21|l in the form: 



ip m (x,uj)(p* m (x',uj) 
2aj a (u - uj m ) 



(29) 



where we have used the fact that uj a ^> uj, u> m . As already 
noted, the {k m = uj m (uj)/c = q m — iK m } are complex but 
for the high-Q modes of interest their imaginary parts, 
K m , will be small compared to the typical spacing be- 
tween their real parts. Therefore the Green function is 
going to be multiply-peaked as a function of uj (which is 
real) around the values u> « cq m . The m term in the 
spectral sum is large when uj = cq m = w m , proportional 
to 1/ n m \ whereas the other terms will have denominators 
at least as large as 5q, the mean spacing between the real 
part of the cigcnfrcquencics (which is roughly the same 
as the mode spacing of the closed cavity). This spacing 
is constant in ID, and only decreases as some power of 
the parameter kR {R is the typical linear size of the res- 
onator) in higher D, whereas the imaginary quantities re m 
can be exponentially small in this parameter. For typical 
dielectric cavity resonators there are many high-Q modes 
for which K m <C Sq and the Green function will be dom- 
inated by the single nearest CF pole at uj = w m . It is 
thus natural to introduce the single-pole approximation 
to the CF Green function in which we replace the full 
Green function near the lasing frequencies by the single 
term involving the nearest pole in the spectral represen- 
tation m. 
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Thus wc assume that the possible lasing frequencies are 
in one-to-one correspondence with the real parts of the 
CF eigenvalues {w m — cq m }. We write uj m w Cl^+Su)^ — 
( K m > 0). In the single-pole approximation, the 
lasing modes are given by ^^(x) = at^ l (pl^ l (x) and are just 
proportional to single CF states. Note that not all CF 
states will lase; the solution Eq. l|2l)|l will find that many 
of the coefficients are zero; the remainder will give the 
intensities of the various lasing modes. Thus, as noted 
above, the lasing modes at this level of approximation are 
identical to (a subset of) the CF states up to an overall 
scale factor. We will work within this approximation 
for the remainder of this article, so we can now adopt 
a simpler notation, dropping the subscript m: — > 

Having just shown that CF states are the correct func- 
tions to describe solutions for high-Q lasing modes of 
open resonators, before we discuss the derivation and so- 
lution of the lasing equations for their intensities and 
frequencies, we will examine some examples of CF states 
and relate them to the more familiar QB states. 



C. Examples of CF states 

1. One- Dimensional CF state 

The simplest example of a set of CF states which also 
has the virtue of being essentially exactly solvable are the 
states of a semi- infinite slab laser (see Fig.QJ. Note that 
this example is also a crude model for an edge-emitting 
semiconductor laser cavity. Consider a dielectric medium 



FIG. 1: Schematic of a dielectric slab cavity of index no and 
length a with perfectly reflecting mirror at the origin. 

of index n — uq which is uniform and infinite in the 
y, z directions and extends from x = to x — a in the 
x direction. At x = it is terminated in a perfectly 
reflecting mirror and from x — a to oo there is vacuum 
(n = 1). To clarify the difference between the CF states 
and the usual cold cavity resonances, first let us consider 
the quasi-bound states of such a resonator. These are 
defined as the solutions of 



-£%tp m (x) = n\\k 2 m Lp m (x) 



(30) 



and the boundary condition (which can moved to the 
interface) d x ip m (x = a) = ik m tp m (x = a). The complex 
eigenfrequencies of this problem are the solution of 



tan(nofc m a) = — in® 
which can be explicitly solved to yield 

n k m a = 7r(m + 1/2) - i- ln[ n ° + ] e 

2 no — 1 



(31) 

n (q m - iK )a, 
(32) 



where m = 1, 2, . . . , oo. Note that the imaginary part of 
the wavevector is always negative, as it should be, and 
in this case is constant, corresponding to the fixed trans- 
missivity of a dielectric interface of index no at normal 
incidence. Since these resonances are outgoing at infin- 
ity, each QB state varies as e +%kmX oc e +K ° x and grows 
exponentially at infinity as mentioned in the general dis- 
cussion of QB states in the introduction. 

Now consider the CF states; inside the medium they 
satisfy exactly the same equation and the continuity con- 
ditions at the dielectric interface are the same as well, 
but the CF state outside the dielectric satisfies the wave 
equation with a fixed external wavevector: 



-dlip m (x) = k 2 ip m (x), 



(33) 



and with the corresponding outgoing wave boundary con- 
dition, d x (p m (x = a) = ikip m (x = a). This defines a 
family of basis states depending on the choice of the real 
external frequency u> = ck. These eigenvalues satisfy the 
equation: 



tan(nofc r „a) = —ino 



k,. 



(34) 



In the Appendix it is shown that this is a well-defined 
eigenvalue problem for any value of k and the eigenvalues 
are always complex with negative imaginary part, similar 
to the usual QB states. But again, unlike the QB states, 
since k is real, these states simply oscillate at infinity as 
e zkx anc j carrv a constant outgoing flux. 

One can see that if k is chosen to be equal to q m , the 
real part of the m th QB state wavevector, then the two 
eigenvalue conditions are almost the same in this vicin- 
ity. Note however that for every choice of the external 
wavevector k there is an infinite set of CF eigenvalues. 
It is shown in the Appendix that from that infinite set 
there is one m th CF state with a wavevector very close to 
that of the corresponding QB state for k = q rn . Specif- 
ically, the difference between the wavevector of the m th 
QB state and the corresponding CF state is given by 

i ~ \ n ° K ° ( \ a<\ 

n [q m - q m )a = — — ; n [K m - K )a = , fe , 2 , (35) 

where /(no) is a number of order unity depending on 
the index, n . Hence this difference is small in the semi- 
classical limit, ka 3> 1. Moreover both the CF and QB 
states are of the form Asin(k m x) within the dielectric, so 
if their complex wavevectors are close, then the CF and 
QB states are almost identical in the medium. Thus we 
see that the standard cold-cavity resonances are almost 
equal to the lasing mode, but not exactly equal. Because 
the CF states are sines of complex wavevector they will 
oscillate and increase in amplitude within the medium 
(Fig.©. 

As noted in the general discussion, these CF states are 
not orthogonal when integrated within the medium. It is 
therefore useful to define adjoint CF states which satisfy 
the complex conjugate differential equation, and with the 
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(a) 3 



(b) 3 r 




FIG. 2: Solutions of the QB and CF state eigenvalue equa- 
tions for a dielectric slab cavity with no = 2. (a) A QB 
state at k m a = 11.51917306 - 0.53647930i. (b) A CF state at 
k m a = 11.55535804 — 0.5310080SH with external wavevector 
ka — 11.51917306. Since the eigenvalues are almost the same 
the two solutions are almost identical within the cavity and 
show large amplification, but they differ qualitatively outside 
as the QB intensity grows exponentially while the CF inten- 
sity remains constant as required for a true lasing solution. 



incoming boundary condition at infinity, d x tp m (x — a) — 
—ikif m {x — a). It is easy to show that the eigenvalues k m 
and eigcnfunctions tp m of this problem are the complex 
conjugates of the CF eigenvalues fe m , and eigenfunctions 
(p m . In the Appendix it is confirmed that these functions 
satisfy the orthogonality relation dxip m {x)(p^ l {x) = 
Jq dx sm(nok m x) sm(nok n x) = r] m (Lo)S mn and the nor- 
malization constant r] m (u>) is calculated. The biorthogo- 
nality of these functions follows from the boundary con- 
ditions at x = a. 



2. Cylindrical CF states 

The next case of interest is the case of a circular (2D) 
or cylindrical (3D) dielectric resonator of uniform index 
no and radius R. The CF eigenvalues can be found by 
applying the conditions Eqs. H24|) and 1|25|) for the circle 
(2D) and cylinder (3D, solutions uniform in z-direction). 
Solutions can be labeled by their good angular momen- 
tum index in the z-direction, M, leading to a countably 
infinite sequence of eigenvalues for each value of M (which 
we will label by m). The states take the form 



(m \JM(n k m r)e ±lM * r<R 

vLKr,M = < JM(n^ H + (fer)e±i M^ r>R 



and the CF eigenvalue condition is found to be 

J' M (n k m R)tt+(kR) = l__fc_ 
U^(kR)J M {n k m R) n k m 



(36) 



(37) 



Note that this equation differs from that defining the QB 
states for this problem simply in that the corresponding 
QB eigenvalue condition would have k m appearing in the 
arguments of H~^, the external solutions, as well as in 
the arguments of Jm, the internal solutions. 

As for the slab resonator discussed above, if we choose 
k> = c<7 m , the real part of a specific QB state frequency, 
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FIG. 3: Solutions of the QB and CF state eigenvalue equa- 
tions for a cylindrical dielectric resonator with no = 1-3 and 
angular momentum quantum number M = 40. (a) A QB 
state at k m R = 71.17700357 - 0.74936509i. (b) A CF state 
at k m R = 71.19539070 - 0.74853902i with external wavevec- 
tor kR = 71.17700357. Again the two solutions are almost 
identical within the cavity; outside the QB solution grows 
while the CF solution decays as 1/r to conserve flux in 2D. 
The amplification of the solutions within the cavity is less 
obvious since both solutions have the usual peak at the caus- 
tic at r/a « 0.45; the decay away from this peak is actually 
slower going outwards then it would be for the cold cavity due 
to the amplification. 



then we will find one CF eigenvalue close to the com- 
plex QB state eigenvalue, as long as the QB state itself 
has small imaginary frequency (high Q). A comparison of 
two such corresponding states is given in Fig. [21 Further 
analysis of the 2D case is given in the Appendix. 



3. CF states for general shapes 

The CF boundary conditions for a general dielectric 
body of arbitrary shape are still easily formulated in 
terms of outgoing spherical waves at infinity, since the 
object will appear point-like at arbitrary large distances. 
We simply require that the CF states satisfy the wave 
equation with index no within the medium and that far 
away they have the form (in 2D) 



¥>m(r 



1 



,</>) = /m(0)-pe 



iki- 



(38) 



In order to satisfy this boundary condition it will be nec- 
essary to solve the wave equation within the body by 
some method and continue it sufficiently far outside to a 
sphere or circle enclosing the body at which radius the 
interior solution can be connected to a superposition of 
outgoing waves with the wave-vector k. The situation is 
simplified if the dielectric body has a smooth shape and 
is not too far from spherical, in which case a Raylcigh 
(Bessel) expansion can be implemented within the cav- 
ity and can be matched to outgoing Hankcl functions 
directly on the boundary specified by R(4>): 



\Em a M 3M(nk m r)c lM * 



r < R(tt>), 
r > R((j>) 

(39) 
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The continuity boundary conditions on the boundary 
can be easily cast into a secular equation and the eigen- 
value condition replaced by a singularity condition. The 
method is then identical to that described in reference [| 
for finding the QB states of such an asymmetric resonant 
cavity (ARC) except that the matrix S which appears 
from the matching conditions is slightly different. Two 
examples from implementing this method for CF states 
are shown in Fig. 0] below. 



IV. GENERALIZED MULTIMODE LASING 
EQUATIONS 




A. Non-linear multimode equations 



FIG. 4: Gray scale plots of CF solutions for a dielectric cav- 
ity with a quadrupolar deformation R{<f>) = Rq{1 + e cos 20) 
at e = 0.16 and an index of refraction n = 3.3 at exter- 
nal wavevector kRo — 20.0. (Left) A bouncing ball mode 
at k m Ro = 20.48417472 - 0.09782270i and (Right) a bowtie 
mode at k m R = 19.71417046 - 0.04025714H. Modes of the 
bowtie type were found to lase in the experiments of Ref. . 



We now return to the formulation of the multi-mode 
lasing equations within the single-pole approximation to 
the Green Function. The outgoing lasing field e(x,u) is 
only non-zero at the lasing frequencies; we approximate 
the Green function in the vicinity of w = f2^ by the term 
arising from the nearest pole, i.e. G(x,x'\oj = f2 M ) « 
(c 2 ^(a;,r2 J11 )^*(£c / ,n /1 ))/(2n 2 a; a (Q A( - Equa- 
tion Q2(J[I now takes the form: 



a fl \l- iD 



1 



dx' 



cavity 

r 



^(x^n^ip^x^n^) 
i + E,<?(aOkl 2 M*',a,)F 



(40) 



where Dq is the scaled pumping rate given by 



non-linear multimode problem can be cast in the form 



D 



Do 



(41) 



We are looking for lasing solutions to these equations, 

i.e. solutions where at least one a M ^ and positive. 

Since a M = 0,V M is always a solution for the steady-state where 

equations, once we reach a pump level at which non-zero 

solutions exist, we must check their stability. 

This can only be done by writing time-dependent equa- 
tions for the CF state amplitudes a M (i) (analogous to the 
standard modal expansions) and determining the effect 
of a small deviation from the coexisting solutions with 
= and a,. 0. Previous work based on the standard 
theory finds [35j that whenever a solution with ^ 
exists, the = solutions are unstable to small devi- 
ations, which then flow to the finite amplitude solution 
with the largest number of lasing modes. We assume 
this property for the current work and will explore this 
stability issue in subsequent work. Ruling out the identi- 
cally zero solutions, and defining I v = |a„| 2 , the general 



ij-Sujfj, + iK^j-iflfj, + 7_l) 



(42) 



1 



dx'- 



i pl{x',n ll )ip ll (x',n IJ ,) 



cavity 



i + J E v g(n v )i v \ip v (x',n v )\* : 

(43) 

and we must look for the solution with the maximum 
number of non-zero positive l v . Note that the CF 
states <Pn(x) in the definition of are understood to 
be evaluated for u> = fi^ and normalized to unity at 
this frequency. We will refer to the full non-linear form 
of the self-consistent equation as the non-linear Haken- 
Saucrmann (NLHS) equations; these equations can be 
treated numerically and we will present results for the 
single and two-mode cases below. 
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B. Generalized Haken-Sauermann equations 

To treat the multi-mode case more simply and to make 
contact with the standard Hakcn-Saucrmann treatment 
we now make the near threshold approximation for 
expanding the denominator of the integrand as 1/(1 + 

E„s(^)k^,(^')l 2 ) ~ i - E,g(^)KM^')\ 2 - we 

then arrive at the generalized (linear) Haken-Sauermann 
(GHS) equations for the modal intensities: 



1 



D J± 



Y,9{n„)A^I v . (44) 



where the overlap matrix is given by 



A tlv = \ dx' ^Ax')^^{x')\^{x')^ 



(45) 



In the limit in which yj_ S> Sl^ and the frequency shift 
Slo^ w Owe recover the standard form of the HS equation 
(Eq. J[J above), with the important difference that the 
overlap matrix A^ v is differently defined; it is complex 
and involves the integral of three outgoing CF states and 
one biorthogonal incoming state, whereas the usual HS 
equations A^ lL , involves the squared moduli of two closed 
cavity states and is purely real. This new form allows us 
to determine iteratively the lasing frequency shifts in the 
multi-mode case, whereas the real approximation does 
not. A second difference is that in the usual HS equations 
the cavity widths are put in by hand at this point, as 
the closed cavity states have zero width. In our approach 
the widths appear automatically as the imaginary parts 
of the complex wavevectors of the CF states. Finally we 
note that in this formulation the theory of interacting 
lasing modes can be seen to have some similarity to the 
theory of interacting electrons in quantum dots, as the 
overlap matrix A^ is similar to a mesoscopic two-body 
matrix element. When the resonator involved has com- 
plex (e.g. wave-chaotic) modes these matrix elements will 
fluctuate and depend sensitively on external parameters 
such as shape and pump profile. Concepts from random 
matrix theory and semiclassical quantum mechanics may 
be useful in analyzing these interactions; an example of 
this approach using the usual near threshold HS equa- 
tions is the work of Misirpashaev and Bcenakker [28( . 



V. SOLUTION OF THE GENERALIZED HS 
EQUATIONS 

A. Iterative method 

Both the non-linear form (|20(l and linear form (|44|) of 
the generalized HS equations require a self-consistent it- 
erative solution since the function F^^I^) in its exact and 
linearized form depends on the unknown lasing frequen- 
cies as well as the unknown intensities I v . Recall that 
the quantities Slo^ and which appear explicitly in the 



NLHS and GHS equations are the difference between the 
lasing frequency f2„ and the nearest CF state complex 
frequency u M (f^ t ); thus both of these quantities are de- 
termined by the real lasing frequencies. Therefore an ini- 
tial ansatz for the lasing frequencies is required in order 
to begin the root search to determine the Con- 
ceptually the linear and non-linear cases are solved by 
similar iterative schemes, the only difference being that 
in the linear case (GHS) the root search becomes equiv- 
alent to inverting a linear system. In the linear case the 
solution must observe the constraint that l v > while in 
the non- linear case the positive roots, if they exist, have 
to be determined (as l v = is always a solution). 

At least two approaches are possible. One can start 
with the N quasi-bound state frequencies which have real 
parts within a linewidth j± of ui a and use Re [uj^ b ] = flff 
as the initial iterate for the real lasing frequencies. A 
second approach is to linearize the eigenvalue equation 
around the atomic frequency tu a , assuming that the imag- 
inary parts of the CF frequencies don't change rapidly 
so that one can use the imaginary parts determined by 
solving the CF eigenvalue condition at u> a as a starting 
point for the iteration. Each of these approaches has ad- 
vantages and can lead to a tractable scheme, as we will 
discuss below in the context of specific examples. We will 
describe the first approach here. 

We first separate the real and imaginary parts of 
Eq. 



1 



1 



£>o7j 



[<5w Jli f2 Jli + k m 7 ± ] = Re [Ffj] 



(46) 



-6uj^± + fi,,^ = Do7j_Im [F/j] (47) 



where we have dropped the tilde on Dq. Our ini- 



(0) 



tial guess for the lasing frequencies is SI 
The quantities Slo^ , f2 M , , F^ in Eq. (|46|) are all func- 
tions of the true lasing frequencies S1 M . However we 
will approximate them by their values at the Sl^ ': let 



5u> 



(0) 



Re 



(0) JO) 



F^F^{ni 0) })^F^;E q 



-Im 
then becomes 



^ M (4 0) ) 



«2V 



Re 



if) ({/„})!. (48) 



This equation can be solved for the modal intensities 1^ 
by varying the pump power and looking for positive roots 
for different trial sets of lasing modes. In practice one 
easily finds the threshold for the first mode to lase and 
then increases the pump in steps above this. 

Having determined the intensities, Eq. H46f) yields an 
expression for the lasing frequencies by substituting Sl^ = 



Soj 



(0) 



.(0) 



and solving for 



(X (°) 

7-l(^ 



D Im 



F, 



(0) 



K (0) n (0) 



1 (°) 



(49) 
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These equations yield the updated lasing frequencies O^ 1 ' 
which can be reinserted into the equations to obtain 
higher iterations. We should note that Eqs. |@8J| and J35J| 
are correct to 0(<5Q M ) and 0(50^) respectively, and we 
have neglected terms of 0(duj^/duj) which can be shown 
to be small in the short wavelength limit. Thus Eqs. 1)48(1 
and 1(49(1 are expected to converge rapidly with iteration 
number. 

The equation for the frequency shifts is a generalization 
of well-known results for the single-mode case and the 
closed cavity. First note that for the closed cavity the 
function Fa is real and doesn't enter the imaginary part 
of the HS equations; hence the pump strength drops out. 
In addition, in a standard treatment the real frequency 
shift 5lo^ is assumed to be zero since one begins with 
the cold cavity frequency, so one obtains: 



AO - It 

osv - . jo) 



(50) 



7-L 



This gives the well known result that when the cavity 
width is large compared to the atomic relaxation rate 



one obtains 5Q L 



(0) 



(0) 



i.e. lasing at the 



atomic frequency u) a ; when << y± (the cavity mode 
is much narrower than the atomic linewidth) one obtains 



-(4 0) mM 0) 



a. 



0^ = Re [uf] , i.e. 
lasing at the cold cavity frequency Our generalization 
shows that openness of the cavity introduces the pump 
strength into these equations and implies a change of the 
lasing frequencies with pump strength. This is a pre- 
diction of our theory which can in principle be tested 
exp er imcnt ally. 



B. Single mode solutions 

1. Near-threshold behavior 

The simplest solution of Eqs. 1)48(1 and 1(49(1 is the case 
where the electric field is oscillating at a single frequency 
fl^. Then, we can solve Eqs. ((48(1 and 1(49(1 straightfor- 
wardly for lu and fi^: 



1 - 



a7 



(51) 



With uniform pumping this equation implies that the 
first mode to lasc (and the only one in the single-mode 
approximation) is indeed the mode with the highest Q, 
corresponding to the CF state with the smallest value of 
its imaginary wavevector k m . (Of course, in a realistic 
scenario the laser typically will not remain single-mode 
for the whole range of pump rates Do and other modes 
can begin to oscillate). As noted above, unlike the closed 
cavity single-mode theory, for which the lasing frequency 



is independent of pump strength, in this case we find 



(0) 



7± A> ff (4 0) )Im 



|(0) 



7-1 



.(0) 



(52) 



The first term in the numerator gives the "center of mass" 
formula for the lasing frequency discussed above, but the 
second term is proportional to 1^ and demonstrates that 
there is a frequency pulling or pushing effect even in the 
single-mode case which depends on the pump strength 
and on the imaginary part of the inverse mode volume, 
A~^. The generalization to multi-mode lasing is here 
trivial with the factor Im [A^^] 1^ simply replaced by 
X^ImLV]!, (cf. Eq. ®). 



2. Far from threshold behavior 

To generalize the above results to arbitrary pump rates 
we simply replace in Eqs. 1(48(1 and (149(1 by its exact 
expression 



f dx'- 

J cavity ^ 



ty l + ^i^OPV 



(53) 



where it is understood that the CF states are calculated 



at the iteration frequency f2 M = fi), 
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FIG. 5: Single-mode single-pole solutions for the modal in- 
tensity vs. pump strength Do for a circular dielectric res- 
onator with an index of refraction n = 2.0 and gain-center 
at uj a = 50.0. (a) Single-mode solution for the intensity I 
of the mode with an angular momentum quantum number 
M = 40 lying closest to the gain center. Dashed line: near- 
threshold solution, solid line: exact solution for modal intensi- 
ties I and dot-dashed line: far-from threshold approximation 
discussed in the text (plotted only for / > 120). (b) Lasing 
frequency Q. as a function of the pump strength Dq. Note that 
Q is measured from its value at the lasing threshold given by 
Qo — 49.6668 and the shift is towards the gain center fre- 
quency. 
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In Fig. [5] we show an example of the proposed iteration 
scheme for a dielectric resonator of circular cross-section. 
There are several observations we can make here. First of 
all, the intensity calculated from the near-threshold the- 
ory greatly underestimates the actual intensity at higher 
pump strengths Do- As we shall further discuss below 
in Section fV CI this will have important consequences for 
mode competition in the multimode regime as the non- 
linear thresholds will be substantially different from those 
predicted by the near-threshold approximation. Second, 
as seen in Fig. |SJb) we now get a frequency shift that is 
power-dependent, with a non- linear dependence on the 
lasing intensities. This dependence is stronger the leakier 
the lasing mode. Third, the exact formula leads to an 
approximately linear dependence of the modal intensity 
1^ on Dq. This linear dependence is termed "satura- 
tion" in laser texts because it is much slower than the 
near threshold rise, but it is still a much stronger de- 
pendence than predicted by the near threshold HS the- 
ory. One can obtain a naive approximation for 1^ by 
replacing 1 + g(f2 M )|( / 9 /J (a;')| 2 / Al in the denominator with 
g(Slp)\<p (l (x')\ 2 I tl , then 

. (**A. ) Do (54) 

where 

* M = f dx' e 2 ^"^ (55) 

J cavity 

where 4>^{x') = — iArg[^ M (cc')]. For comparison, we plot 
the result of this expression in Fig.[5jb) (dot-dashed line). 
We see that this formula is qualitatively correct but over- 
estimates the modal intensity. 



not present realistic calculations on complex cavities us- 
ing our generalized formalism; we defer such calculations 
to future work. Instead we present a simple compari- 
son of the near-threshold and non-linear lasing solutions 
for the two-mode case, which confirms that the widely- 
used near threshold approximation can introduce a large 
quantitative error in the lasing thresholds and hence over- 
estimates greatly the number of lasing modes for a given 
pump power in the complex dielectric cavity lasers we are 
studying. We solve Eq. (|42f) for two interacting modes of 
the dielectric cylinder (other modes which might also lase 
are neglected). In Fig. [S]and[7|we show the results. For 
mode 1, with the lower threshold, the near-threshold and 
non-linear calculation give the same threshold as they 
must (Do ?s 0.11), and then deviate substantially well 
above threshold as we already saw in the previous sec- 
tion. Mode 2, which would have a threshold of Dq ps 0.18 
without mode competition has a threshold of Do as 0.25 
in the near threshold approximation, including the effects 
of mode competition, but has a correct threshold, taking 
non-linear effects into account, of Do ~ 0.41. Hence we 
see that in the large interval 0.25 < Dq < 0.41 the near- 
threshold theory predicts that two modes would be lasing 
and with comparable intensity, while the full theory pre- 
dicts only one lasing mode. Moreover, because of this 
delay, when Mode 2 does begin to lase, its intensity lags 
mode 1 by an order of magnitude in contrast to the pre- 
diction of the near threshold analysis. We expect these 
discrepancies to be enhanced if there are many compet- 
ing modes and many fewer modes to be lasing at a given 
pump power than predicted by the near-threshold theory. 
These issues will be explored in detail in future work. 



VI. SUMMARY AND CONCLUSIONS 



C. Two-mode solutions 

The numerical effort required to solve the linear and 
non-linear generalized HS equations for a non-trivial case, 
e.g. a dielectric cavity with fifty modes underneath the 
gain curve, is substantial. Moreover the interesting ques- 
tions require varying the pump strength and the shape 
of the cavity. In earlier work we have done extensive cal- 
culations along these lines using the near threshold HS 
equations for the closed cavity and then using reasonable 
but uncontrolled approximations to evaluate the modal 
intensities outside the cavity [2i|. The current general- 
ized formalism gives us a method to evaluate emission 
directionality and output intensities correctly and with- 
out any approximations beyond the usual ones allowing 
a multi-mode solution. It also presents a formulation of 
the lasing equations valid away from the threshold, using 
the general function, F(I V ), defined in Eq. l|4^|) above. A 
striking shortcoming of the near threshold approximation 
that we found in our earlier work psj is that it appears 
to overestimate greatly the number of lasing modes com- 
pared to the experimental observations [t| . Here we will 



We have reformulated semiclassical lasing theory to 
treat open cavities, with particular attention to complex 
and random lasers, for which the output power and di- 
rectional emission patterns are not trivially found from 
knowledge of the internal modal intensities. We begin 
from the assumption that a steady-state multi-periodic 
lasing solution exists and is stable. A key idea is to for- 
mulate the theory in terms of the self-consistent linear 
response to the polarization described by an outgoing 
Green function. This yields a set of self-consistent inte- 
gral equations for the lasing modes which can be solved 
by various means. The approach which is closest to con- 
ventional multi-mode lasing theory is to write a spectral 
representation of this Green function in terms of outgo- 
ing states with constant flux. For high Q cavities these 
states are similar but distinct from the usual resonances 
or quasi-bound states; they satisfy useful biorthogonal- 
ity relations with adjoint functions. The Green func- 
tion is then approximated by a single-pole near each las- 
ing frequency implying the lasing mode is just an in- 
tensity 1^ times a single CF state. General and near- 
threshold equations then follow for these intensities 1^ 
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FIG. 6: A single-pole two-mode solution for the modal in- 
tensities as a function of pump strength Do for a circular 
dielectric resonator with an index of refraction n = 2.0 and 
gain-center at uj a = 20.0. We have chosen 7x = 50 using 
effectively a flat gain curve. The two modes have the angu- 
lar momentum quantum numbers Mi — 20 and (b) Ms = 18 
respectively. The modal intensities Ii and I2 are given in (a) 
and (b). Solid lines denotes the exact solution, while we plot 
the near threshold approximate solution in dashed lines. The 
non-interacting thresholds are given by = 0.1069 and 

Dth = 0.1767 (marked with vertical arrow in (b)). By def- 
inition modal interactions cannot change the first threshold, 
which is thus the same in the near-threshold and exact so- 
lution. However interactions dramatically modify subsequent 

thresholds; in this case the second threshold is increased to 

(2) 

Dll' = 0.2464 in the near threshold approximation, but ac- 
tually should be increased much more, to Dl h ' = 0.4139, 
based on the exact solution. This indicates that even taking 
modal interactions into account through the near-threshold 
approximation is inadequate and likely greatly overestimates 
the number of lasing modes at a given pump power. 

and for the lasing frequencies Q^. The near threshold 
equations are generalizations of the well known Hakcn- 
Saucrmann equations derived for a closed cavity. These 
equations can account for the effects of mode competi- 
tion and spatial hole-burning which are important but 
difficult to treat in complex cavity lasers. An iterative 
scheme was described for solving these equations and 
some illustrative results were given for the single-mode 
and two-mode case. The one and two-mode solutions 
indicate that the near-threshold approximation substan- 
tially overestimates the number lasing modes well above 
threshold and underestimates their output power. 

In our view this approach clarifies the longstanding 
question of how to treat rigorously open cavities within 
semiclassical lasing theory. It now remains to apply this 
formalism to cases of interest: dielectric micro-cavity 
lasers of interesting shape, cavities with fully chaotic ray 
motion, and lasing from random media. 



APPENDIX A: PROPERTIES OF THE 
CONSTANT FLUX BASIS 

In this Appendix we will explore the constant flux (CF) 
states, as defined in Sec. IIII C II We investigate the an- 
alytically tractable cases of the ID dielectric slab laser 
and the 2D cylindrical laser, and argue that the qualita- 




FIG. 7: The corresponding laser frequency variation with 
pump power for the example in Fig. |^|as a function of pump 
strength Do, for the exact solution. The frequencies of the 
two modes are measured with respect to their values at their 
respective lasing thresholds given by f2i — 19.8584 (full line) 
and 0,2 — 20.3764 (dashed line). Note that since the gain 
center is uj a — 20.0 both modes are pulled towards the gain 
center. 

tive features of the solutions for these simple geometries 
should also hold for more general shapes. Here we fo- 
cus on three properties of the CF basis. First, that for 
each high-Q QB state of the cavity there exists a single 
CF state with similar real and imaginary part and hence 
very similar behavior within the cavity (where they sat- 
isfy the same differential equation with just a slightly 
shifted eigenvalue). Second, we prove that all the CF 
states have eigenvalues with negative imaginary parts, 
corresponding to amplification of the wave within the 
cavity. Third, we confirm explicitly the biorthogonality 
of the two (outgoing and incoming) CF bases within the 
cavity for a general geometry. 

1. Dielectric slab cavity 

As stated in Sec. 1111 C II for a semi-infinite dielectric 
"slab" resonator the QB and CF eigenvalue equations 
are: 

tan (nk rn a) = —in (Al) 

and 

tan (nk m a) = —in ^ — - (A2) 

Where k m = q m — ik m and k m = q m — in rn are the 
complex QB and CF wavenumbers, a is the length of 
the slab, n is the (possibly complex) index of refraction, 
and k is the wavenumber corresponding to the external' 
Fourier frequency. As in the main text we use tildes to 
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denote quantities associated with the QB states, but we 
drop the subscript on the index of refraction. We also 
define the dimensionlcss quantities x m = nq m a, y m = 
and z = nka. Using trigonometric identities, 
we can expand both equations into real and imaginary 
parts: 



and 



cotq m (l — tanh 2 k, 
cot 2 q m + tanh 2 k r , 



tanh K m (cot q m + 1) 
cot q m + tanh K m 



cot x m (l — tanh 2 y„ 
cot 2 x m + tanh 2 y„ 



(A3) 



(A4) 



(A5) 



the spectral representation of Eq. (|29|l . We are therefore 
especially interested in finding solutions to the CF eigen- 
value equation with real part close to the wavevector, 
k, appearing as a parameter in the CF eigenvalue equa- 
tion. Furthermore, since we expect the real part of the 
QB mode frequencies to provide a first approximation to 
the actual lasing frequencies, we look for solutions to the 
CF eigenvalue equations with the external Fourier fre- 
quency set equal to the real part of a quasibound mode 
frequency. Equations (|A3() - (|A6() suggest that, for each 
QB mode with a sufficiently short wavelength, a corre- 
sponding CF solution exists with approximately the same 
eigenvalue, since in this limit x = nk m a is large. In that 
case the right hand side of (|A5|) is small so that (|A5|) is 
nearly identical to l|A3jl . and for q m w x equation (|A6(I 
approaches i|A4(l . 



tanh y m (cot 2 x m + 1) 
cot 2 x m + tanh 2 y m 



(A6) 



In the method used to solve the multimode laser equa- 
tions presented in this paper, the Green's function is ap- 
proximated by a single pole, i.e. by a single CF state in 



The solution to Eqs. 1)A3|) and l)A4j) was given in the 
text by equation l|32|) . and equations (|A5|) and (| A6|) can 

be solved for tanh y m and cot x rn : 



tanhy m = 



-(cot 2 x m + 1) ± \J (cot 2 x m + I) 2 - 4(nx m /z) 2 cot 2 



(A7) 



CO tj X fii — 



(1 - tanh 2 y m ) ± J (1 - tanh 2 y m ) 2 - A(ny m /z) 2 tanh 2 

2ny m /z 



(A8) 



Note that in order to self-consistently choose solutions to 
(|A7|) and (|A8|) which are close to QB solutions, we must 
select the minus signs in both equations. We now expand 
15|) in the small difference 5q m = q c f — q q y. 



6x m {l-l/n 2 ) 
1/n 2 



Assuming that 6y m is, at most, first order in the small 
parameter z , we can replace y m with k m , leading to: 



Sx m (n — 1) s=s k 



(A10) 



Or, using the expression for ko from Eq. (|32|) and the 
definitions of x m and z: 



Sq n 



L n — 1 J 



i 



(All) 



2(n 2 — 1) na 2 k 

We see the difference between real part of the CF reso- 
nance and the external frequency is small when (1) the 



external frequency is set near the real part of a QB res- 
onance and (2) the semiclassical parameter z _1 is small. 
To find an approximation for K m , as well as to justify the 
approximation used in deriving i|All|) that y m w k m , we 
next expand (| AT|) to second order in Sq m x and x -1 . After 



(A9) some manipulation we obtain: 



Sk Ti 



2(n 2 - l) 2 

f(n) 1 
(ka) 2 a 



inl^i] 1 
n — 1 \ 



In[2±i] 

In— 1 J 



i 



na 3 k 2 
(A12) 



Which is the result quoted in Section IIII C II Thus we 
see that Sq m is first order and 6K m is second order in z^ 1 . 
Fig. JHJl shows that the approximation used here agrees 
well with the numerically evaluated results for the differ- 
ence between the QB and the nearest CF wavevector. 
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(a) 




60 80 100 120 140 160 



(b) 




40 60 80 100 120 140 160 

q m 

FIG. 8: The difference between the CF and QB eigenvalues 
for a ID dielectric slab resonator, as a function of the real part 
of the QB eigenvalue, (a) The fractional difference between 
the real part of the eigenvalues, and (b) the fractional differ- 
ence between the complex parts. In both plots, the triangles 
represent the actual difference, found numerically, which is 
compared to the linear approximation, Eq. I|A120 . derived in 
the text and represented by the solid curves. 



2. Cylindrical dielectric cavity 

It is well known that for a cylindrical cavity, there exist 
whispering gallery mode solutions to the QB eigenvalue 
equation that have very small widths, k m . Moreover, it 
is these modes which will dominate lasing emission from 
dielectric cylinder lasers. In this section we show that 
there exist corresponding solutions to the CF eigenvalue 
equations close to these QB solutions. Specifically, the 
difference between the CF and QB lifetimes is propor- 
tional to a small factor times the QB lifetimes; hence 
these CF and QB states have almost identical lifetimes. 

The QB eigenvalue equation for a cylindrical resonator 



r(+), 



^-J M {nk m R)l 'J M (nk m R) = ^-H { +\k rn R)/ H [ A J> (k m R) 
ar dr 

(A13) 

Within the regime of interest where the modes are well 
trapped (and the index of refraction is not too small), 
we can approximate the Hankel and Bessel functions ap- 
pearing in this equation by the Debye approximations 

1371. Define tanh a = 



{nk m R/Mf 



1. 



l-(fc m i?/M) 2 and tan/3 = 
In terms of these variables, the De- 



bye approximations are: 

JM(k m R) = 

Y M {k m R) = 



g M(tanh a-a) 

\J2-kM tanh a 



- Af(tanh a— a) 

i-TrMtanha; 



J M (nk m R) 



Y M (nk m R) 



cos (M tan/3 - MP - f) 



\-kM tanh a 



sin (M tan /3- Mp- \ ) 



r7rMtanh< 



(A14) 



(A15) 



(A16) 



(A17) 



Note that in the regime in which these approximations 
are valid, we can ignore the r-dependence outside the 
exponent when we evaluate its derivatives, so that the 
eigenvalue equation is approximately: 



tan P tan A ~ tanh a 



2ie 



(A18) 



and 



where we have defined A = M(tan/3 — /?)—§ 
H = A/(tanha; — a). By expanding in the small quantity 
k m and neglecting e M compared to e _/i (which is valid 
for well trapped modes), we can derive an approximate 
expression for k m : 



M tanh a 
n 2 a 2 q m (n 2 - I)' 



(A19) 



A similar, but not identical, result was obtained by 
Nockcl j23j , where it was pointed out that the small factor 
e 2fJ - could be interpreted as the exponential decay factor 
due to tunneling. 

The eigenvalue equation for CF modes in a cylindri- 
cal resonator has the same form as Eq. (|A13ll above 
for QB states except that the functions a, p, A and 
fi appearing in the CF equation are evaluated at dif- 
ferent arguments. Explicitly, A and P are evaluated at 
k m R = k m R + S(k m R), whereas \i and a are evaluated at 
the external Fourier frequency, which we have set equal to 
q m . We can then expand in the small quantities d(k m R) 
and k rn . To first order, the imaginary part of this equa- 
tion yields: 



(A20) 



M{n 2 - 1) tanh a' 

where it is understood that tanh a is evaluated at 
tanh a = 



1 — (^gp) 2 Note that in the semiclassical 
regime whispering gallery modes have a high M value, 
so the difference between the lifetimes for the CF and 
QB mode is a small factor times the already small QB 
lifetime k„, . 
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(a) 



In Kn 



(b) 





FIG. 9: (a) The log of the imaginary part of the QB resonance 
eigenvalue for a cylindrical resonator with no — 2, plotted 
against the real part of the eigenvalue; triangles represent 
numerical results compared to the analytic approximation of 
Eq. 1A19L (b) The fractional difference between the CF and 
QB lifetimes versus the real part of the QB eigenvalue for 
the same resonator. Stars are numerical results, solid line the 
approximation of Eq. IA20II . We have plotted the sequence 
corresponding to angular momentum M = 50. 



3. General argument that Im(k m ) — 



< 



We were able to show above explicitly that for the CF 
eigenvalue closest to the corresponding QB state eigen- 
value, K m > for both the ID slab laser and the 2D 
cylindrical laser. Here we show that this result holds gen- 
erally for all CF eigenvalues and for arbitrary geometries. 
This is important because the exact CF Green function 
involves all CF eigenvalues and it would be unphysical 
for it to have poles in the positive imaginary plane. 

The equation satisfied by the CF states is: 



nPh^cpfa) (inside) 
k 2 if(x) (outside) 



(A21) 



We define a vector flux as / = ^-(i^* Vv? — ipVip*). The 
outgoing wave boundary conditions, as expressed in the 
main text by equation (|24|l . imply that if we integrate the 
radial component of / over increasingly larger spheres S 
(or circles in the 2D case) with origin within the cavity, 
the value of this integral approaches a positive constant 



proportional to k: 



fdA 



(A22) 



Note that V ■ / vanishes outside the cavity, whereas inside 



it is equal to —ri (p(p*Im(k^ l ) 



2n 2 \ip\ 2 q r , 



By using 



the divergence theorem and the CF equation we can also 
write this integral as: 



f dA= V • fdV = 2n 2 q m K m / n 2 \(p\ 2 dx 

J Int{S) J cavity 

(A23) 

Comparing equations l|A22(l and (|A23|) we can see that in 
order to satisfy l|A22|) . which follows from the outgoing 
wave boundary condition, we must have K m > 0. 



4. Biorthogonality of incoming and outgoing CF 
states 

Consider a solution to the CF eigenvalue equation ip m 
with eigenvalue fc m , along with an adjoint solution Tp n 
with eigenvalue fc*. Using the fact that these two func- 
tions satisfy the eigenvalue equations (|22|l and l|26|l . re- 
spectively, we have: 

n o( k l - fc m) / dxlp* n ip m 

J cavity 



/ dx {ip* n V 2 tf m - (^ m V 2 (,o*) 

J cavity 
J boundary 



(A24) 



Where the last integral is over the surface of the cavity 
boundary. Outside the cavity we have V • (Tp^Vip m — 
ipmVtp^) = 0, so we can in fact replace the integral over 
the cavity boundary with an integral over a large sphere 
S inclosing the cavity. We can then invoke the outgo- 
ing/incoming wave boundary conditions to evaluate the 
surface integral explicitly: 

/ (K^™ - VmV^*) • dA 
Js 

— > / {i-WkVm - ikVm<?n) ' dA = (A25) 
Js 



For the case k n ^ k m , equations (|A24|) and l|A25(l there- 
fore imply: 



d xip* m ip n = 



(A26) 



cavity 



It should be noted that, although the CF functions are 
defined to exist both inside and outside the cavity, the 
inner product is defined by an integral over the cavity in- 
terior only. The biorthogonality of the CF states is there- 
fore determined completely by their behavior within the 
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cavity and at its boundary; we merely use the behavior 
of the CF functions at infinity to show that the surface 
term in equation l)A24(l vanishes whenever the CF states 
satisfy the outgoing/incoming boundary conditions. 

We can check the biorthogonality of the CF states ex- 
plicitly for the case of the dielectric slab cavity, and also 
calculate the normalization factor ?y m (u>) defined in equa- 
tion l|28|) . In the cavity interior, we have 

ifn(x) = sin(nofc„a) (A27) 



and 



ip* m = sm(n k m a) (A28) 



Using the trigonometric identity sin(s) sin(t) = i(cos(s — 
t) — cos(s + 1)), we can write the inner product (for the 
case m ^ n) as: 



sin(no k n a) sm(nok m a)dx 

cos(nofc ra a) cos(nak m a) I tan(nofc n a) — tan(nofc m a) tan(nofc„a) + tan(nofc m a) 
2tiq y k n k m k n + k m 



(A29) 
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